Size-dependent fine-structure splitting in self-organized InAs/GaAs quantum dots 
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A systematic variation of the exciton fine-structure splitting with quantum dot size in single 
InAs/GaAs quantum dots grown by metal-organic chemical vapor deposition is observed. The 
splitting increases from -80 to as much as 520 /ieV with quantum dot size. A change of sign is 
reported for small quantum dots. Model calculations within the framework of eight-band k p theory 
and the configuration interaction method were performed. Different sources for the fine-structure 
splitting are discussed, and piezoelectricity is pinpointed as the only effect reproducing the observed 
trend. 



The exchange interaction of electron-hole pairs (exci- 
tons) in semiconductor quantum dots (QDs) has been 
subject of a lively debate in recent years 0, 0, S S IE 
lI M ■ In such strongly confined systems it is supposed 
to be enhanced with respect to the bulk case due to the 
close proximity of electrons and holes. However, the in- 
fluence of the exact geometry of the confining potential 
on the exchange interaction still needs to be clarified. A 
detailed understanding of the resulting exciton fine struc- 
ture in quantum dots is of fundamental interest and of 
largest importance for potential applications of QDs in 
single-photon emitters and entangled two-photon sources 
for quantum cryptography 

The total angular momentum M of heavy-hole excitons 
(X) in QDs is composed of the electron spin (s = ±|) 
and the heavy hole angular momentum (j = ±§), con- 
sequently producing four degenerate exciton states fre- 
quently denoted as dark (M — ±2) and bright (M = ±1) 
states indicating whether they couple to the photon field 
or not. Independent of the given confinement symmetry 
electron- hole exchange interaction causes a dark-bright 
splitting. Furthermore it mixes the dark states lifting 
their degeneracy and forming a dark doublet (|2)±|— 2)). 
Likewise, additional lowering of the confinement symme- 
try to C^ v or lower mixes the bright states producing a 
nondegenerate bright doublet (|1) ± | — 1)). 

While emission lines involving pure states are circularly 
polarized, the mixed states usually produce lines show- 
ing linear polarization along the [110] and [110] crystal 
directions, respectively (Fig. QJ. The two bright states 
are thus directly observable as linearly polarized transi- 
tions in luminescence experiments. The energetic differ- 
ence between these lines is called exciton fine-structure 
splitting (FSS). 

The biexciton (XX) ground state is not split by the 
exchange interaction, since the net spin of the involved 
electrons and holes is 0. However, the XX to X decay in- 
volves two allowed transitions with the final states being 
the bright states of the X. Therefore, the FSS is repro- 
duced (yet inverted) in the XX to X decay (Fig. ^| . 

Recently reported experimental values of the FSS in 
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FIG. 1: Energy scheme illustrating the effect of the exchange 
interaction on the exciton (X) bright states. XX denotes the 
biexciton ground state. \M) is the exciton total angular mo- 
mentum, a+/— denotes circularly polarized light, and ir x / y 
linearly polarized light with x = [110] and y = [110]. D2d, 
Ci v , and C2v indicate the confinement potential symmetry. 



In(Ga)As/(Al)GaAs QDs are listed in TablclJ Its magni- 
tude has been shown to be influenced by annealing 5, 6], 
and external electric fields 0. A systematic analysis of 
the influence of the geometry of the confining potential 
and an explanation for the wide range of observed ener- 
gies however, is still missing. 

Although most of the cited experimental references 
in Table provide some theoretical background, only 
Bester et al. Q linked the FSS to morphological QD 
properties in a theoretical work. Using an empirical 
pseudopotential-based approach including the configura- 
tion interaction (CI) they analysed In^Gai-^As dots of 
various shapes and composition. The largest FSS they 
obtained was 30 /iteV for an elongated QD with a lat- 
eral aspect ratio of 1.3. They pointed out that even a 
cirular based QD having a structural Cm symmetry ex- 
hibits a non-zero FSS due to the atomistic asymmetry of 
the underlying lattice. In their work, however, they did 
not account for the piezoelectricity-induced asymmetry, 
which plays a major role in our interpretation of the large 
FSS observed. 

Major effects that lead, in principle, to a lowering of 
the confinement symmetry include (a) structural elonga- 
tion of the QDs, (b) strain-induced piezoelectric fields 
|10|. and (c) interfacial symmetry lowering and its en- 
hancement by atomistic elasticity 0, |n| . Since the lat- 
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TABLE I: Examples of experimental values of the exciton fine- 
structure splitting (FSS) in quantum dots. Column 1 gives the 
material system and column 3 the experimental setups used 
(fjPh = microphotoluminescence, FWM = four- wave mixing) . 
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ter effect has been shown to produce FSS values smaller 
than 10 fxeV we concentrate on the two former effects 
in our modeling. 

We present in this Letter a systematic study of the FSS 
in single InAs/GaAs QDs by experiment and theory. 

The sample investigated was grown by metal-organic 
chemical vapor deposition on GaAs(OOl) substrate. A 
300 nm thick GaAs buffer layer followed by a 60 nm 
Alo.6Gao.4As diffusion barrier and 90 nm GaAs were 
grown. For the QD layer nominally 1.9 monolayers of 
InAs were deposited followed by a 5 s growth interrup- 
tion. Subsequently, the QDs were capped with 50 nm of 
GaAs. Finally, a 20 nm Alo.33Gao.67As diffusion barrier 
and a 10 nm GaAs capping layer were deposited. Anti- 
mony was added during the deposition of the InAs layer 
and the subsequent growth interruption. 

The sample was examined with a JEOL JSM 840 scan- 
ning electron microscope equipped with a cathodolumi- 
nescence setup ^3 ■ It was mounted onto a He flow cryo- 
stat which provided temperatures as low as 6 K. The 
luminescence was dispersed by a 0.3 m monochromator 
equipped with 1200 or 1000 lines/mm gratings. Above 
1.2 eV the light was detected with a liquid-nitrogen 
cooled Si charge-coupled-device camera and below 1.2 eV 
with a liquid-nitrogen cooled InGaAs diode array. At 
1.28 eV the minimal spectral resolution as given by the 
setup was «140 fieV. 

In order to reduce the number of simultaneously de- 
tected QDs a metal shadow mask with circular apertures 
of 100 and 200 nm in diameter was applied onto the sam- 
ple surface. This technique allows the reproducible inves- 
tigation of single QDs. 

The QD ensemble luminescence peak is centered at 
1.18 eV with a FWHM of « 150 meV. It displays a 
pronounced modulation which is due to a shell-like size 
distribution of the QDs producing a number of different 
subensembles. The QDs display a large substate splitting 
suggesting In/Ga interdiffusion to be negligible [Lj- In 
addition, transmission electron microscopy (TEM) im- 
ages taken from similar samples indicate a truncated 
pyramidal shape |14| . Therefore the exciton energy is 
directly connected to QD size. 

When measured through the apertures the lumines- 
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FIG. 2: Polarized spectra for two different QDs emitting 
at high and low energies are shown. The FSS is —80 and 
420 /ieV, respectively. Gray lines in the lower panel origi- 
nate from charged excitonic complexes not considered in this 
Letter. 

cence decomposes into single sharp lines. We are able 
to unambiguously identify spectra of single QDs consist- 
ing of up to six lines via their characteristic spectral dif- 
fusion pattern 15). From these lines the exciton (X) 
and biexciton (XX) transitions can be extracted by per- 
forming excitation- and polarization-dependent measure- 
ments. Thus the FSS in the corresponding QD can be 
directly measured. Figure [21 shows examples of polarized 
X and XX emissions for two different QDs. 

The FSS shown in Figgis -80 ficV for the QD emit- 
ting at high energies and 422 /j,eV for the QD emitting 
at low energies. The FSS is defined to be positive if the 
X line at lower energy is polarized along the [110] crystal 
direction. The order of the X and XX lines is reversed 
due to a varying impact of correlation in the QDs as has 
been established earlier [lfij . 

Figure shows measured FSS values as a function of 
the exciton recombination energy, i.e., the center between 
the two X lines. Over the whole energy range of QD 
luminescence a systematic variation of the FSS within the 
very same sample is observed. Small QDs that emit at 
high energy show the smallest FSS, while large QDs that 
emit at lower energies show the largest FSS. This trend is 
surprising, since the dark-bright splitting has been shown 
to become larger the smaller the QDs are p| . 

The largest obtained value amounts to 520 /xeV, 
exceeding all previously reported FSS values in the 
In(Ga)As/GaAs material system, where the FSS is re- 
ported to be always smaller than 180 /ieV (Tab.P). Only 
Finley et al. observed in the InAs / AlGaAs system values 
up to 1 meV 7] . They proposed the stronger confinement 
due to the AlGaAs matrix as a possible reason. 

As the FSS is related to a lateral anisotropy of the con- 
fining potential, we investigate the two most promising 
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FIG. 3: Measured exciton fine-structure splitting (FSS) as a 
function of exciton recombination energy for 49 QDs. The 
gray curve depicts the QD ensemble luminescence with num- 
bers denoting QD height in units of InAs monolayers. A clear 
trend of the FSS is observed. The data points extracted from 
Fig- HI are indicated by red (gray) circles. 



candidates leading to such an effect, namely an elonga- 
tion of the QD along [110] or [110] direction and strain- 
induced local piezoelectric fields. A thorough theoreti- 
cal investigation requires two essential ingredients: First, 
precise information about the QD structures involved 
and, second, a theory that reliably links the QD structure 
with their electronic and optical properties. 

Structural assumptions: In Refs. |13|, LLJ] it is demon- 
strated that the QDs are of pure InAs and have the shape 
of truncated pyramids with abrupt interfaces. The mod- 
ulation of the low-intensity-photoluminescence (PL) (see 
Fig. |2l directly allows to determine the exact height of 
the QDs belonging to different subensembles. We con- 
sider both square-based QDs and, although plan-view 
TEM data show no indication for a major QD elonga- 
tion [TJ(, QDs elongated along the [110] crystal axis. The 
base lengths of the model QDs are chosen in a way that 
the exciton spectrum matches our PL and PL-excitation 
spectra ^3| and TEM data The smallest square- 

based QD has a base length of 10.2 nm and a height of 
3 monolayers, the largest one a base length of 15.8 nm 
and a height of 13 monolayers. Elongated QDs were de- 
signed to have the same volume as their square-based 
counterpart with the same height. 

Method of calculation: The excitonic states are cal- 
culated using the CI method by expanding the exciton 
Hamiltonian into a basis of antisymmetrized products of 
single particle wave functions (6 electron and 6 hole wave 
functions), thus accounting for direct Coulomb interac- 
tion, correlation, and exchange |l5j. The single particle 
orbitals are derived from a strain-dependent eight-band 
kp Hamiltonian accounting for band coupling and band 
mixing, following our previous work |16|. Absorption is 
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FIG. 4: The piezoelectric potential calculated for (a) the InAs 
bulk and (b) the extrapolated ei4 value. The potential is 
shown for a lateral cut near the base of the QD (white square) , 
which is modeled as a truncated pyramid having a base length 
of 12 nm. 



calculated by Fermi's golden rule applied to CI states. 

To assess the role of piezoelectricity quantitatively, an 
accurate value of the prime input parameter, the piezo- 
electric constant en must be known. However, as Bcstcr 
et al. pointed out recently > the standard reference for 
the InAs piezoelectric constant eu = 0.045 C/m 2 0] 
cannot be considered a trustworthy value for strained 
InAs QDs. More recent measurements [l8l fl9l I20I l2l| 
on strained InGaAs quantum wells grown on GaAs(lll) 
substrate suggest an ei4 value with a different sign and a 
much larger magnitude. Because of the apparent strain 
dependence of ei4, the exact value for strained InAs 
QDs is unknown. To encounter this problem we extrap- 
olated the values of the most recent work |2l| to ob- 
tain the value of InAs. Remarkably, the obtained value 
ei4 = —0.385 C/m 2 has a different sign than the 'stan- 
dard' value for bulk InAs as has been noted before by 
Bester et al. ^lj. Since there is no direct experimental 
confirmation for this value, quantitative results of the cal- 
culations have to be treated with caution (see below). In 
addition, we point out that the exchange interaction was 
calculated including only monopole terms. Calculations 
including dipole terms are under way. 

The impact of the 'new' en value compared to the 
standard value on the piezoelectric potential is threefold 
(see Fig. 0J: First, and of most importance for the dis- 
cussion in this Letter, the magnitude of the piezoelectric 
potential inside the QD is strongly enhanced by a fac- 
tor of around ten. Second, the field orientation inside 
the QD is rotated by 90° and third, the charge distri- 
bution near the QD corners change from dipole-type to 
monopol-type. 

Although both piezoelectricity and elongation intro- 
duce an anisotropy of the confining potential, their im- 
pact on the orientation of the participating electron and 
hole wave functions is fundamentally different: Piezo- 
electricity forces electron and hole groundstate wave 
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FIG. 5: Calculated exciton fine-structure splitting (FSS) val- 
ues for a series of QDs with a square base including piezoelec- 
tricity (bl) and a series of QDs elongated in the [110] crystal 
direction in absence of piezoelectricity (b2). C(n,n) denotes 
the number of electron and hole states included in the CI cal- 
culation. The structural assumptions of the modeled series 
are shown in (al) and (a2). 

functions to point in orthogonal directions, whereas the 
groundstate wave functions located in an elongated QD 
accommodate their orientation to that of the QD, thus 
being parallel to each other. 

To distinguish the impact of elongation and piezoelec- 
tricity, we calculated the FSS for a series of square-based 
QDs including the effect of piezoelectricity [Fig. [5] (al)] 
and a complementary scries of QDs elongated along [110] 
with a lateral aspect ratio of 2 in absence of the piezo- 
electric field [Fig.[5](a2)]. 

Fig. [5] (b) shows the resulting FSS values for both se- 
ries. The obtained trend for series 1 (piezoelectricity) is 
in qualitative accord with the experiment (Fig.|3J- The 
importance of the piezoelectric potential for the discus- 
sion of the FSS stems from its strong dependence on 
QD height, which has already been shown earlier [l6j . 
independent of the chosen value of ei4. This explains 
why the FSS strongly increases with increasing QD size 
(decreasing exciton energy) [Fig. Efbl)]. Quantitative 
agreement, however, could not be reached and deserves 
further investigation. A strong elongation of the QDs in 
the [110] direction by contrast does not reproduce the ex- 



perimentally observed trend [Fig. GJb2)]. An elongation 
in the perpendicular [110] direction (which is equivalent 
to reflecting the results on the x axis) also fails. It gives 
strongly negative FSS values for small QDs and slightly 
positive values for large QDs. Therefore elongation can 
be ruled out as the source of the observed FSS trend in 
the QDs investigated. 

For small QDs the influence of piezoelectricity is weak 
and the FSS scatters around /xeV (Fig. |3J). Then, sec- 
ond order effects like small random elongations of the 
QDs or interfacial symmetry lowering govern the FSS. In 
some QDs this leads to negative FSS values. 

A considerable scatter of the observed FSS values of 
more than 100 /ieV for a given energy is present in Fig. [3] 
It is due to variations in QD geometry for a given QD 
height. A different base length for example leads to a 
different strain distribution inside the structure which in 
turn leads to different piezoelectric fields and hence to 
different FSS values. The variation in QD geometry for 
one QD height is visible through the observed FWHM 
of the subensemble peaks of around 30 meV in the PL 
spectrum 0]. Another source for the FSS scatter can 
be slight in-plane elongations which are more or less pro- 
nounced. 

To asses the influence of correlation on the FSS we 
compared calculations with a CI basis of 6 electron and 
6 hole levels (C(6,6)) thus including correlation effects 
to calculations with only 2 electron and 2 hole levels 
(C(2,2)) excluding correlation. Qualitatively the results 
are the same, eliminating a decisive influence of correla- 
tion on our conclusions (Fig. [5j|. 

In conclusion, for the first time a systematic change of 
the FSS with QD size is observed. We found a variation 
by one order of magnitude and observed a sign inver- 
sion of the FSS. Large values of up to 520 /ieV for large 
QDs were found. Complementing model calculations us- 
ing eight-band k-p theory and CI rule out elongation of 
the QDs as a source for the large FSS. To assess the influ- 
ence of interfacial asymmetry in detail, atomistic theories 
need to be considered. Piezoelectricity can explain the 
size dependence but quantitatively fails to explain the 
huge variations observed. 
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